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Abstract 

We introduce a method to obtain the analytic solution of the higher-order Baxter equation for 
twist-two and twist-three operators of planar A/" = 4 SYM. Our result proofs the conjectured 
formula for the three-loop anomalous dimension of twist-two operators. As such we derive the 
maximally transcendental part of the corresponding three-loop QCD result from the maximal 
supersymmetric gauge theory in four dimension purely by methods of integrability. 
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1 Introduction 



The exploration of the AdS / CFT correspondence using the methods of integrabihty has led 
to a first insight into its interpolating property between weak and strong coupling. Since the 
discovery of integrable structures on both sides of the correspondence [T] many techniques 
have been introduced and developed which will hopefully lead to the complete solution of the 
planar spectral problem in a finite volume. 

Integrable structures in four-dimensional quantum field theories are not new and have 
also been found in non-supersymmetric theories before, see for example [21 [3l H] and [5] and 
references therein. However, superconformal invariance preserves integrabihty of the maximal 
supersymmetric planar A/" = 4 SYM to higher, possibly all, orders. The factorized scattering 
of elementary excitations of a super-spin chain [6j describing = 4 SYM at one-loop is 
extended to higher-loop order [7J and allows the determination of the asymptotic spectrum 
in terms of long-range Bethe ansatz equations [H] • The ambiguity of a phase factor of the S- 
matrix, determined by symmetry solely [9], is also fixed by a constraining crossing-symmetry 
[lOj supplemented by an unequivocally specification of its solution pT| [T^. 

The operators which gained the most attention so far are twist operators with lowest 
possible anomalous dimension, whose thermodynamical behavior is governed by the so-called 
universal scaling function [121 [T3] . It was shown that the strong coupling expansion of this 
scaling function, see ^4] and [15] and references therein, correctly reproduces the know string 
results [16] of the corresponding dual state to two-loop order in perturbation theory. 

To this class of operators belong the shortest possible local composite operators of the 
A/" = 4 theory, the twist-two operators. For a simple representative of these one starts from 
the protected half-BPS states Tr Z"^ and inserts M covariant derivatives T) 

Tv [ZV^^ Z) + ... . (1.1) 

In the spin chain picture this is a non-compact s[(2) spin = — ^ length-two Heisenberg 
magnet with M magnons. The dots indicate the mixing of all states where the covariant 
derivatives may act on any of the two fields. For each even M there is precisely one highest 
weight non-BPS state whose total scaling dimension is 

oo 

A = 2 + M + 7(5,M), with ^{g,M)=Y,l2i{M)g^' , (1-2) 

1=1 

where 7(5, M) is the anomalous part of the dimension that depends on the coupling constant 
= and A = N g^^^ is the 't Hooft coupling constant. 

For twist-two closed expressions for the anomalous part 7(17, M) of the dimension are 
known to two-loop order from explicit field-theory calculations [17j, and at three-loops from 
a solid conjecture [18] extracted by the principle of maximum transcendentality [19] from the 
QCD splitting functions at three-loops [ID]. To this order the anomalous part 'y{g, M) of the 
dimension can also be reliably computed by the asymptotic Bethe ansatz [7] for fixed values 
of M. Up to relatively high values of M it was checked ^ that it coincides with the two- 
and three-loop anomalous dimensions of the twist-two operators, which are known in terms 
of nested harmonic sums as obtained in [TTJ [18] . Due to wrapping effects that appear at four- 
loop order, these operators have also been used to show the incorrectness of the asymptotic 
Bethe ansatz at this order of perturbation theory [21j . 
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It had not been known how to obtain these closed results in a rigorous analytical way from 
the Bethe ansatz. This problem will finally be solved in the present paper. By assuming the 
principle of maximum transcendentality it is possible to make a general ansatz that contains 
harmonic sums of the proper degree with unknown coefficients and to fit these with a high 
numerical precision from the Bethe ansatz, but this principle is not granted in general for 
other natural operators [22] . 

Therefore, we introduce a systematic approach to derive these results from the Baxter 
equation that can be obtained from the s((2) Bethe ansatz. As such we derive the maximally 
transcendental part of the corresponding QCD computation from the maximal supersymmet- 
ric Yang-Mills theory in four dimensions solely by means of integrability. As an application 
of our method we also derive the three-loop anomalous dimension of twist-three operators. 

The method presented in this paper is based on deforming the one-loop solution. The 
precise structure of this deformations can either be obtained by taking the Mellin transforma- 
tion of the perturbative Baxter equation or simply by an educated guess, based on comparing 
the higher-loop Baxter equation with the difference equation of a general deformation of the 
one-loop solution. In the main text we follow the latter, while the aspects of the solution in 
Mellin space are given in the appendix. 

We begin our analysis by introducing the perturbative Baxter equation, which has been 
proposed in [2l|. After constructing the solution for twist-two operators in an pedagogical 
way, we will apply the same techniques to twist-three operators but only state the final result. 
Lengthy and detailed computations are shifted to the Appendix as well as the complete 
calculation of anomalous dimension. We close with a summary and give some ideas of further 
applications. 



2 Baxter equation 

The Bethe roots that parametrize a solution of the Bethe ansatz equation obey the following 
expansion in the coupling constant up to a loop order i 

uu{9') = Y. g'' uf . (2.1) 

i=0 

They are given by zeros of the corresponding Baxter function Q{u) which also exhibits a 
similar expansion 

M £-1 

Q{u) = n - = E Q^'^ ■ (2-2) 

k=l i=0 

To three- loop-order Q{u) satisfies the Baxter equation in form of a second-order finite differ- 
ence equation [21] 

A+(a;+) Q{u + i) + A„(x~) Q{u - i) = tL{x) Q{u) . (2.3) 
The variables are defined through the spectral parameter u in the following way [25] 



x-(n) = - 1 + . (2.4) 
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To the first three orders of perturbation theory, t = 2, the functions Aj- have the following 
form 

A±(x) = exp ( + ^7£^ " 4:^ ^ogQ^'\n) \ . (2.5) 



^ X ^ X ^ dn^ 
The auxiliary transfer matrix tL{x) for twist-two and three operators is given bj0 

tL{x) = 2u^ + q2{g^)u^-^ , i = 2,3, (2.6) 

where the charge q2 is also to be expanded in g^. Using (j2.2p one can read off the coefficient 



functions , see below. The resulting expressions for the anomalous dimensions are given 

by 

7(9') = 7+(5')-7-(5') (2.7) 

and 

£-1 



2 i=0 

3 Twist-two operators 

We will now explain how to find the higher order Baxter function for operators that have been 
analyzed to a large extend in AdS/CFT, twist-two operators. After restating the one-loop 
problem with its known solution which for completeness is also given in Appendix El we focus 
on the explicit construction of the higher loop solution. The main idea of the construction 
is based on deforming the one-loop solution in a suitable manner. As mentioned before, one 
can also perform the complete computation in Mellin space and verify the result, as included 
in Appendix 1X1 

The leading order Baxter equation is given by 

(u + §)' g(°)(^/ + i) + (n - §)' QW(n - i) - (u)QW(n) = , (3.1) 
with the leading order transfer matrix given by 

if {u) = 2u^ + ) , = -M{M + 1) - i . (3.2) 

The solution to this equation is given in terms of the continuous Hahn polynomial 



1 . (3.3) 



g(0)(,) = 3^^(-M,M + l, i + ^n 

It is rather straightforward to compute the anomalous dimension using the above explicit 
formula for the Baxter function 

7(°)=8 5i(M). (3.4) 

At next-to-leading order the Baxter equation can be obtained by expanding the corresponding 
asymptotic all-loop equation, see [21], and is given by 

(u + §)' Q«(n + i) + {u- §)' Q^^\u -i)- tf{u)Q^^\u) = (3.5) 
= (2-i ^-^{u + I)) gW(n + i) + (2 + i 2j^(^x - §)) QW(n - i) + t^^\u)Q^''\u) , 



^Note, that ga is zero for the twist-three ground state ^l] 
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with t 



(1) 



4 + ^■y^^\2M + 1)). In order to solve the above equation one is lead to consider 



the following class of continuous Hahn polynomials 

y{u) = 3F2 



n, n + a + b + c + d— 1, a + iu 
a + c, a + d 



(3.6) 



which satisfy the difference equation |26) 

n{n + a + h + c + d-l) y{u) = B{u) y{u + i) - [B{u) + D{u)\y{u) + D{u) y{u - i) (3.7) 
where 

B{u) = (c — iu){d — iu) , D[u) = (a + iu){b + iu) . 



By a deformation of the one-loop solution (jS.Sp we understand a suitable choice of a, 5, c and 
d, such that when the deformation parameters are set to zero, the function (j3.6p coincides 
with ()3.3p . One of such deformations is given by 



-M, M + 1 + 25, \ + iu 



1 + 6,1 

Upon differentiating once w.r.t. 5 aX, 5 = Q the resulting difference equation reads 

[u + Q^'{u + i) + {u-^f Q^'iu -i)- tf{u)Q^'{u) = 
= -i{u + §) Q^^\u + i) + f(n - §) Q(°)(n - f) - (2M + 1) 1 

In the above formula we have used the following notation 
d 



(3.8) 



)(0) 



(n). 



86 



Qf{n) 



Q^'{u) and Qf{u) 



<5=0 



<5=0 



Thus using the linearity of (j3.5p one can multiply (|3.8p with to obtain the part of Q^^\u) 
which is proportional to the one-loop anomalous dimension. 

The missing parts of the full solution can be found by considering the following deformation 
of the one-loop solution 



Q 



(0) 



(3.9) 



(3.10) 



3-f2 



-M, M + 1, i + ZM 



(3.11) 



1 + 6, 1-6 

The second derivative of the resulting difference equation evaluated at 6 = takes the form 

{u + §)' Q^"(n + i) + {u-^f Qf{u -i)- tPiu)Q^"{u) = 2Q^^\u) - 2Q(o)(n + i) . 

(3.12) 

To complete the solution one must note that a third type of the deformation must be intro- 
duced such that the argument of the r.h.s. of the above equation is shifted u — s- (u — i). From 
(|3.7p one infers that this is obtained by interchanging a ^ c and b ^ d. Thus, the last term 
reads 



3F2 



-M, M + 1, ^ + m + 6 



1 + 6, 1 + 6 



(3.13) 



Since the one-loop solution (|3.3p is also a solution to the homogeneous part of the two-loop 
Baxter equation (13. 5p one can add it to Q^^\u) with an arbitrary coefficient function of M. To 
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fix the function uniquely one notices that if the leading order Baxter function is a polynomial 
of degree M, then the higher loop corrections are polynomials of degree M — 1 and for even 
distribution of roots, as in the case considerecH, M — 2. Hence, one has to add the term 
a(M) Q^^\u), where —a{M) is the ratio of the highest order term u*^ of Q^^\u) and Q^'^\u) 
given by 

a{M) = A{S2{M) +ASi{Mf - 2 5i(M)5i(2M)) . (3.14) 

Finally, the full next-to-leading order Baxter function with the appropriate normalizatiorH i 
given by 



IS 



g«(n)=a(M)3F2 



-M, M + 1, i + m 



2 86^ ^ 



86^ 



86^ 



1, 1 



-M, M + l + 2d, 1 + iu 
1 + 6,1 

-M, M + 1, l + iu 
1 + 6, 1-6 

-M, M + 1, '^+iu + 6 
1 + 6, 1 + 6 



5=0 
1 



5=0 



(5=0 



(3.15) 



With these results it is rather straightforward to compute the two-loop anomalous dimension 
in a closed form, see Appendix |Bj It is given by the following known expression, where all 
harmonic sums are evaluated at argument M 



7 



(1) 



(M) = -16 (S3 + S-3-2 S_2,i + 2 Si{S2 + 5-2)) . 



(3.16) 



At three-loop order one finds the following perturbative Baxter equation for Q^'^\u) 

{u + ^)^Q^^\u + i) + {u-l)^ g(2) _ _ tp (u) g(2) (u) = 

= (2-i^{u + §)) Q«(n + i) + (2 + i^{u- I)) Q«(tx - i) + 4'^(tx)Q«(u) 
+Piu) iu + i)+ P* (n) (n - i) 
-{^i2M + 1) + 2^2 + {u + i)"' + (n - i)"')Q(°)(^z) , (3.17) 

where P{u) and K2 = K2{M) are given by 

(3.18) 

It is instructive to split the solution into five different sub-classes, which are independent 
of each other due to the transcendental composition of the corresponding inhomogeneities 
and can be analyzed separately as presented in Table [TJ 

■^This follows from expanding Uk{g^) in (12. 2p in a power series in g. 

^The global normalization of the full Baxter function Q{u) = Q'''^\u) Q^^\u) can be chosen arbitrarily. 
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Index 


Classification 


I 

n„ 
lib 

lie 

III 


polynomial parts of Qo 
terms of order )0 
terms of order (7'''^-')^ 
terms of order (7*-*^'*)^ 
non-polynomial parts of Qo 


IV 


one-loop normalization 



Table 1: Classification of solutions 



The classes I, IIq, II^ and lie can be obtained by using the same arguments as for the 
next-to-leading order analysis. Terms polynomial in Qq lead to 



7W+a(M)7W d 
2 85 
_K2{M) + 2a{M) 

2 W 

K2{M) + 2a{M) d"^ 
2 5^ 

The Ila part of the solution is given by 



3-^2 

3-^2 



-M, M + 1 + 25, \ + iu 
1 + 5, \ 

-M, M + 1, \ + iu 
1 + 5, 1-5 

-M, M + 1, i + m + 5 
1 + 5, 1 + 5 



<5=0 



5=0 
1 



5=0 



^6 9(54^^ 



92 92 ^ 

+ 77-0 -7777^ 3^2 



M, M+1, i+m 
1 + 5, 1-5 

-M, M + 1, i + m + 5 
1 + 5, 1 + 5 

-M, M + 1, i + m + /3 
l + a + P, l-a + P 



5=0 
1 



Polynomial inhomogeneities linear in 7^°^ redound to the result 



5=0 
1 



a,/3=0 



= — 



(0) 



^3 



The term of ©((7^'^^)^) is given by 

^3lIcN-^— ^3i^2 



-M, M + l + 2a, ^+ra + /3 
1 + « + /?, 1 + /3 

-M, M + 1, ^ +ZU + /3 
1 + /3, 1 + /3 



a,/3=0 



-M, M + 1 + 2Q;, i + iu 
l + a + /3, 1-/3 



-M, M + 1 + 25, ^ + i« 
1 + 5, 1 



/3=0 

1 



a,/3=0 , 



a=o 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



To find the solution for class III, these are all the terms that contain the denominators 
{u±i/2) ' needs to proceed more generally. The difference equation satisfied by continuous 
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Hahn polynomials does not allow to generate (^^"^^^2) terms by the deformations previously 
introduced. However, all terms of type III add up to a polynomial and one can make use 
of a general statement formulated in lemma 1 to find the solution, see Appendix O where 
also details of the calculation can be found. When brought together, the terms in class III 
can be expressed by a sum of two polynomials. The first part is given by a hyper geometric 
polynomial, and the corresponding solution can be easily found using the techniques discussed 
in this section 



Q 



(2) 

Ilia 



U 



7 



(0) Q2 Q 



4 da'^dp 

^(0) Q2 Q 



4 da^ 8/3 



-M, M + 1, i + m, 1 
1 + a, l-a + P, 1-/3 

-M, M + 1, \-iu, 1 
1 + a, l-a + P, 1- P 



Q:,/3=0 



a,/3=0 



(3.23) 



On the other hand, the second polynomial is much more complicated and to find the solution 
one has to apply lemma 1 from Appendix ICl 



Qliil(^) 



E 

fc=0 



mkR{k,M) 



{{^ +iu)k + - iu)k) , mfc = 



aj~i (j - 1)! 

1 fRij,M) ■ 



(3.24) 



See App. [A] and App. [C] for the definitions of R{j, M) and aj. Finally the solution to the 
homogeneous equation, i.e. the leading order solution with the correct normalization is fixed 
again by imposing the proper degree reduction of Q2{u) such that Q2{u) is a polynomial of 
degree M — 2. The normalization is computed in the same fashion as in the two-loop case. 
Hence is given by 



Q\t^iu)=a^^\M),F2 



-M, M + 1, i + 



1, 1 



(3.25) 



with 



a(2)(M) = -(525i(M)^ - 5_3(M)(95i(M) - 65i(2M)) - A8Si{MfSi{2M) 

-35_2(M)(75i(M)2 - 45i(M)5i(2M) + S2{M)) + 3Si{Mf{ASi{2Mf 
-bS2{M) + AS2{2M)) - 6S4(M) + ?,{S2{Mf + 25i(2Af)(53(M) - 25_2,i(M))) 



-25i(M)(ll53(M) -95_2,i(M))) 



2m M ■ 



(3.26) 



The anomalous dimension obtained from this result in terms of nested harmonic sums with 
argument M is given by 

7(2) (M) = 64(25_5 + 2S5 - 45_4,i - 25_3,-2 - ^-3,2 - 25_2,-3 - 85i,_4 - 4Si,4 

-95'2,-3 - 5^2,3 - 253-2 - 553,2 + 2( - 254,1 + 5'_2,-2,l + 5'_2,l,-2 + 45i,_3,i 
+ 'S'l,-2,-2 + 'S'1,-2,2 + 65i,i,_3 + 25i,i,3 + 25i,2,~2 + 25i,2,2 + 25i,3,i + 352,-2,1 



+2(52,1,-2 + 'S'2,1,2 + 'S'2,2,1 + 'S'3,1,1 — 25i,i,_2,i))^ • 



(3.27) 



It agrees precisely with the conjecture of [18], upon changing the basis of the harmonic sums. 
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4 Twist-three operators 

A special class of maximal helicity twist-three operators in QCD has also found to be inte- 
grable to the first loop correction [27j and the relation of the solution to Wilson polynomials 
has been established by the Baxter approach in [28j. 

In this section we will apply the methods of the previous one to obtain the three-loop 
expression for the Baxter function of twist-three operators of = 4 SYM 



Tt {VZV^ ZV'^ Z) + ... , si + S2 + ss = M 
The leading order Baxter function is given by the Wilson polynomial [2 



gW(n) 



4-^3 



_M M,^ 1 , ■ 1 



m 



1, 1, 1 

from which one computes the anomalous dimension 



(0) 
73 



8Si(f). 



The Baxter equation for the twist-three operators ()4.ip to two-loops is given by 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



+ 3 



(0) 



+ ) (n + §)QW(n + i) + (3 + i2^{u-^)) [u - ^) Q^'^^u - i) , 



(0) 



where 



2^3 - (M^ + 2M + |)u, 4 (^) = (-73 + (4.5) 
Using the same techniques as for the twist-two case one finds the final result to be given by 

1 I 1 



Q«(^.) = c(M)4F3 



' 2 ' 2 + ^' 2 ~'~ 2 



m 



1, 1, 1 



(0) 



3 92 



1-^3 



f , f + 1 + 25, i + ^ 
1 + 6, 1 + 6, 1 

2 ' 2 ^ ' 2 ^ ' 2 



5=0 



2 86"^^ ^ \ 1 + 6, 1-6, 1 

with the normalization function c{M) given by 



(5=0 



c(M) 



„(o) 



T3 



25i(M)) +352( 



(4.6) 



(4.7) 



The anomalous dimension computed from these closed expressions, see App. IB. 21 is given by 



7l')(M) = -8( 53(f) + 25i(f)52(f 



(4. 



as has been guessed in [2T1 [29] . 
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The NNLO correction to the Baxter function can be found as the solution of 



(4.9) 



v(o) 



-{u 



+ 1)^ (^i + i)QW(n + «) + (^3 + ^ 



„{o) 



-(n 



u 



+ 4'^(^)Q^'n^^) + ^3(^^) Q^^'Hu + + P^{u) g(°)(n - i) - {-fPiM + 1) + 2K3)g(°)(tx) 
where -P3(u) and = K^[M) are given by 



+ 1)' ,K3(M) 



(4.10) 



Following the same method outlined in section [3] and considering each part of the solution 
separately according to table [T] it is straightforward to find the complete solution. The 
normalization function c^^^ for the homogeneous solution is again fixed by degree reduction 
i.e. Q(2) 

is a polynomial of degree (M — 2), 



c(2)(M) = 325i(f )^ - 32Si{ffSi{M) + ^S2{ff + 8Siiff{Si{Mf + S2{M)) 
+4Si{M)S,{f ) - 45i(f )(5i(M)S2(f ) + 353(f)) - |S4(f ; 

The three-loop Baxter function is given by 

.M M + i 1 + i 
1, 1, 1 

M M 



(4.11) 



Q(2)(^)=c(2)(M)4F3 



m 



7^ ^ + c(M)7f a 
+ 4 9^^^^ 

/^3(M)+3c(M) 



1 + 5, 1 + 5, 1 



2 ' 2 ^ ' 2 ^ ' 2 



3 ^ 



A/ M 

1 + 6, 1-6, 1 



2 ' 2 ^ ' 2 ^ ' ' 2 



l + (5, l-(5, 1 

iu 



5=0 



5=0 



3 (0) d 



1 



(0) 



+ 



+ 



9a a/32 

92 

9a2 5^ 



1 + a, 1 + Q + /3, 1-/3 
iu 



Q,/3=0 



32 952 



(7: 



(0)n2 



92 



32 952 



2 ' 2 ' 2 ' 2 
1 + a, l-a + /3, 1-/3 

■f, f + 1 + 25, i + in, i 
1 + 5, 1 + 5, 1 



in 



5=0 



2 ' 2 ^ ' 2 ^ 2 
1 + 5, 1-5, 1 



lU 



(4.12) 



5=0 



For the three- loop anomalous dimension one obtains, again see Appendix IB. 21 for details, 

7(2) (iW) = 8(55-25i,4-652,3-10S3,2-654,l + 8(5i,2,2 + 52,l,2 + 52,2,l + Sl,3,l + S3,l,l)), (4.13) 

with all sums evaluated at f . The result coincide^ with the conjecture of [2T| [29]. 



Note, that we choose to write the harmonic sums in the canonical basis. 
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It is important to note, that the three-loop Baxter equation for the twist-three operators, 
in contradistinction to the case of twist-two operators, does not contain superficially non- 
polynomial parts and thus the solution can be found in terms of deformations of the one- 
loop solution only. Although the wrapping problem for twist-two operators starts from the 
four-loop order, this superficial breakdown of the polynomiality of the Baxter equation may 
signalize its incorrectness at the next order. The same applies for the twist-three operators, 
where the corresponding four-loop Baxter equation contains rational functions, though the 
solution should correctly reproduce the corresponding anomalous dimension [21j. 

5 Summary and Outlook 

We have shown how to solve the two- and three-loop Baxter equation for a special subset 
of operators. With the explicit form of the Baxter function we were able to reproduce the 
known results based on Feynman calculus [TT] and to prove the three- loop conjecture [18] 
for twist two-operators that has e.g. been used to check the field theory solution for the 
three-loop planar dilatation generator obtained by algebraic methods in [30]. Likewise, we 
gave a proof for the anomalous dimensions of twist-three operators that were conjectured in 
\21\ [29] . Besides our approach to find analytic solutions to the Baxter equation there are 
also techniques to obtain such solutions directly from Bethe ansatz equations, see [22] and 
especially [31] and references therein. 

Due to the nature of the mechanism one can trace back all different contributions to 
the anomalous dimension to the corresponding inhomogeneity of the perturbative Baxter 
equation. As such it would be of great importance to generalize the successfull application 
of TBA [3^ of the four-loop Konishi [33], i.e. twist-two M = 2, operator to all twist-two 
operators. The knowledge of the correct four-loop result of twist-two operators in a finite 
volume could than be used to analyze and, if possible, to fix the four-loop Baxter and hence 
Bethe equations. A detailed application of our methods to the next order should also reveal 
a different structure appearing in the Baxter equation that renders the twist-two solution 
incorrect but, in turn, should give the right result for twist-three. 
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A Solution in Mellin space 

In this appendix we will show how to derive the one and two-loop solution of the Baxter equa- 
tion for twist two operators in Mellin space following [2]. We consider the Mellin transform 
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Q^'^\lj) of the Baxter function Q^^\u) 

POO 

Q^)[u) = / (iL^c^™-iQ»(cj). (A.l) 

In what fohows we will make use of the relations 

Q(^) (u + ai) = doj u}""'^ \uj-^Q^^ (w) } , 

/■CXD / d \ J C 1 

{u + biYQ^'\u + ai)= duJ''-^Lo-^{iuo—^^J-''Q^^{u)Y (A.2) 



One-loop. The differential equation we obtain from the leading order Baxter equation (13. ip 
reads 

By means of the variable transformation z = 1/(1 — w) and by applying the Faddeev- 
Korchemsky substitution, see [4], Q(*)(u;) Q^'^\—oj) in (|A.ip we replace (|A.ip and (|A.3p 

by 

— ■ /"^ 



Q^-^(n)= . / dz(l-zr-^z— ^QW(^), (A.4) 

r(iM)r(i - iu) Jq 

respectively. It is convenient to introduce new functions Q^'^\z) 

Q»(z) = 7i(r^gW(z), (A.6) 
for which ()A.5j) takes a simpler form 

)(A)%a-24-,f-i}a'»'(.)^0, (A.7) 

where we have omitted the factor ^J z{l — z). It should be noted, that both Q(^)(z) and 
q(o)(1 _ satisfy the above equation. Being a differential equation of the second order, 
equation (|A.7p has two algebraically independent solutions and both, Q^'^\z) and Q^'^\l — z), 
can be written as linear combination of these two basis vectors. In particular, one may 
impose symmetry properties on the solutions. We choose the two independent solutions to 
obejH Q(0)(i _ ^ Q(0)(z) and Q^^\l - z) = -Q%\z). The first solution can be found with 
the ansatz 

oo 

q(o)(^) = ^ (A.8) 

fc=0 

with arbitrary coefficients Ck and a. The beginning of the resulting recurrence leads to the 
consistency condition a = and comparing terms with same powers of z, we find the relation 

{k + lfCk+i = {k-M){k + M + l)Ck. (A.9) 



^In what follows wo assume even values of M. 
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This first order recurrence is solved by 

Ck = Rik,M)Co^ 

where 



R{k,M) 



T{k + M + l)T{k - M) {-l)T{k + M + 1) 



(fc!)2r(M + l)r(-M) (A;!)2r(M + l-A;) 
Thus, one solution to the Baxter equation ()A.7p reads 



g(o)(^) = ^0 2F1 



-M, M + 1 
1 



(A.IO) 
(A.ll) 

(A.12) 



Using ()A.4p we obtairlf] the following final result in the variable 



r(m)r(l - iu) 



-M, M + 1, i + m 
1, 1 



(A.13) 



Note, that the factor containing gamma functions is just a phase and hence negligible. 

The second solution can be found by including logarithmic terms in addition to (jA.Sp . 
Imposing the antisymmetry, we are looking for a solution of the following form 



<5b (z) = Yl (""k log z - at log(l -z)+hk) R{k, M) . 



(A.14) 



A:=0 



Terms of power log(z) and — log(l — z) lead to the relation 

(A; + if R{k + 1, M) ak+i = {k - M){k + M + I) R{k, M) , (A.15) 
which by definition of R{k, M) in (jA.lip simply leads to 

a-k+i = ak = ao. (A.16) 
All terms of z^ lead to the intertwining relation for the coefficients and bk 

2ik+l) R{k + 1, M) ak+i + {k+lf R{k + 1, M) b^+i = {k-M){k+M+l) R{k, M) bk , (A.17) 
which can be solved straightforwardly 

bk = bo-2aoSi{k). (A. 18) 

Please note, that one can make use of the antisymmetry of the solution to fix bo to be 

bo = bo{M) = 2Si{M)ao . (A.19) 
Upon integration of ()A.14p by means of taking a derivative of the Beta integral 



Jo k\ 



(A.20) 



^For the integration we make use of the invariance Q^'^\z) = Q'"'(l — z 



1(0)^ 
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the second solution is given by 



M 



Qs i^) = 2^ {l)k{l)kk\ ^ +iu + k)- ^0(2 



fe=0 



tu) 



2Si{k) + 2Si{M)) . 

(A.21) 

The relation of this solution to the one already obtained can be seen using the following 
identity for the Polygamma function 



It can then be written in the following form 



d 



-M, M + 1, l + iu + 6 
I + 6, 1 + 6 



lU) 



in tanh(7rn). 



(A.22) 



<5=0 



+ (25i(M) +mtanh(7ru)) Q'^^^u). 



(A.23) 

The first part of this solution agrees with (jA.lSp when multiplied by the normalization factor 
—2Si{M) while the second part is nothing but (|A.13p multiplied by a periodic function and 
a coefficient function of M. 

Note, that there is still another representation for Q^'^\ Making a different choice of 
variables p = 4z(l — z) to rewrite (|A.5p and performing the same steps of the computation 
one finds 



T{iu)T{l-iu) 



1> 1, 



(A.24) 



This representation is equal to (IA.13p . see for example 

Two-loops. The differential equation that corresponds to the two-loop Baxter equation (13.5 
for the same choice of variables ()A.6p reads 
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(0) 



d 



- — (1 - 2z)— + 2M] + — 

2 ^dz I z{l-z) 



(A.25) 



As the acting differential operator is linear we will analyze the two inhomogeneous terms 
separately, starting with the part proportional to 7^^^^. It can be rewritten as 



where the coefficients -B^ are given by 

{k - M){l + k- M) 



By 



k + l 



R{k, M) Co 



(A.26) 



(A.27) 



To find the first part of the solution Q^j^\z) we make the ansatz 



= E^fc^(^'^)■ 



(A.28) 



k=0 
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which leads to the condition 



7^°) 1 + k-M 



The solution is given by 



Dk = Do + 



7 



(0) 



[2 Si{k + M) - 2 Si{M) - Si{k)) Cq . 



(A.29) 



(A.30) 



Thus, the first part of the result can be written in a compact form, noting that 



d (1 + M + 25)a 



dS (1 + 6)k 
it is given by 



<5=o 



(2 Si{k + M)-2 Si (M) - Si (k)) 



(A.31) 



Q^')(z) = Z)o 2i^i 



-M, M + 1 
1 



<5=0 



(A.32) 



In n-space the result is obtained according to (jA.4p and hence reads 



Q'-j\u)=ADo 3F2 



-M, M + 1, l + iu 



1, 1 



-M, M + 1 + 25, l + iu 
1 + 6,1 



(A.33) 



5=0 



(A.34) 



with, for completeness, the phase A given by 

^ i7rr(i + iu)r(i — iu) 
r(m)r(l - iu) 

Finally, let us focus on the second term to complete the solution, i.e. 
|z(l - z) + (1 - 2z)^ + M(M + 1)| g«(z) = -^^^ E ^) • (^-35) 

One realizes that for this specific differential operator, we have to make the ansatz that 
~ ^og^ z Y^^Lq z'^ to obtain any expression that contains 1/z terms. However, the 
resulting expression also includes terms of order log^z, logz, z'^ such that we need three 
consistency conditions to meet the requirement of the r.h.s. of (|A.35p . Therefore we should 
consider an ansatz of the following form 
00 

Qb^ (^) = E ("'^ (log ^ - log(l - ^))' + (l°g ^ - - ^)) + ^'^) ^) • (^-36) 

fc=0 

Comparing the powers of log^ z, log z and respectively leads to the following recurrences 
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" fc + 1 
2 

Cfe+i = Cfe - ^fc+i and ao = -Cq 



(A.37) 
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The solutions to ()A.37p are given by 



bk = bo -AaoSiik), Ck = cq - 2bo Si{k) + 8ao 5'i,i(A;) . 



(A.38) 



Note, that (|A.36p is the natural transcendental generalization of the second solution of the 
leading order Baxter function (IA.14|) and as such its 6o(l ~ (/c)) part leads to (|A.23p and 
is not of importance since the one-loop part will be fixed in the end by requirement of degree 
reduction. For the same reason it is also dispensable to fix cq by symmetry requirements. In 
addition writting the double index sum as 2Si^i{k) = Si{k)'^ + S2{k) the result is given by 



M 



(-M)fe(l + M)fc(i + iu) 



A;=0 



{l)k{l)kk\ 



+ iu + k) + ^i{\- iu) + 2 S2{k) 



+ (*o(^ + iu + k)- ^'o(i - in) - 2 Si{k)f + 2 Siik)) . 



(A.39) 



One of the terms containing 2S2{k) can be written in a compact form by absorbing it into a 
deformation of the Pochhammer symbol 



1 



d5^ {l + 5)k{l-5), 



5=0 



(i)fc(i); 



■S2{k)^ 



(A.40) 



Likewise as in the case of the second solution at leading order we use the identity ()A.22h and 



m) 



to obtain the following form 



cosh (iru) 



-M, M + 1, l + iu + 6 
1 + 6, 1 + 6 

-M, M + 1, i + m 
1 + 6, 1-6 



5=0 



5=0 



(A.41) 



(A.42) 



We have neglected all one-loop parts and phases and fix the overall leading order influence to 
the two-loop Baxter function by the fact that its degree should be M — 2. The final result is 
given in (|3.15p . 



B Anomalous dimensions 

To compute the anomalous dimension from the closed expressions obtained, it is important to 
note that the twist-two and twist-three Baxter functions are real and hence invariant under 
the map u — > —u. Furthermore the leading-order Q-functions evaluated at +^ are equal to 
one, i.e. Q^^\+^) = 1. 

B.l Twist-two operators 

Two-loops. From ()2.8p one infers the following form of 7^^^ 

j(^\m) = A{M) + B{M) , (B.l) 
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where 

A{M) = i( + 2 QW'(n)3 - 3 Q(°)'(^/) Q(°)"(n)) "^""^ = ^+ - A_ , (B.2) 

S(M) = 2if Q(^)'(n) - Q(i)(n) ^("^'(u)') = B+ - B_ . (B.3) 

It is convenient to express Q^^^ (u) in the fohowing way 

Qi^\u) = J:'-^C, + ^u),, with h,M = i-lf('!)(''^'). (B.4) 



Taking the derivative and using some identities for nested harmonic numbers, see [35], one 
finds the closed expressior|j| for the terms entering A_ (M) 



M 



Q{o)'(_i) = i 5^^^^^ Si{k) = 2i 5i(M) , (B.5) 

k=l 
M 

=i2 ^ 6^^^ (2^1,1 (A:) - 2^2 (A:)) = 4^2(2^1,1 - ^2 + S.2) , (B.6) 



k=l 
M 



g(0)///(_ I) ^ -3 ^ ^^^^ (6Si,i,i(A;) - 65i,2(A;) - 6^2,1 (A;) + 653 (A;)) 



k=l 



24 (2 5i,i,i - Si,2 - ^2,1 + 5i,-2 - 5-2,i) • (B.7) 



Hence, j4_ is given by 



A_{M) = 8(35_3(M) - SsiM) - 65_2,i(M)) . (B.8) 
One of the sums that needs to be evaluated for i?_ is 

Q«(-i) = a(M) +4 5_2(M) -8 5i,i(M) +4 52(M). (B.9) 
To find the second part we are in need of the following auxiliary formula 

M 

J2 bkM Siik + M) Si{k) = (-1)^'(8 5i,i(M) - 5 52(M)) . (B.IO) 

k=l 

After some algebra Q^^^'l-f) is given by, harmonic sums evaluated with argument M, 

Q(i)'(_i) = 2i{ a{M) 5i - 4 (2 53 + 5_2,i - 5i,_2 + 5i( 5-2 + S2) 

-3(5i,2 + 52,i-2 5i,i,i))), (B.ll) 

such that reads 

S_(M) = 16( 53 - 5^3 + 5i( 52 + 5_2) + 2 5_2,i) • (B.12) 



^Please, be reminded that all states of s[(2) have even M. 
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Due to the symmetry of the Baxter function, = —A- and = — -B-, the two-loop 
anomalous dimension of twist- two operators is given by (j3.16p . 

Three-loops. It is straightforward to expand the formula for 7^^^ in terms of the perturbative 
Baxter functions. In order to find its close form we need to evaluate the following sums of 
the leading-order result 

M 

Qm^)^_^) = iiY,KM (24Si,i,i,i(A:) - l2Si{kfS2{k) - l2Si{k)) 

k=l 

= 96 i^(2(2S'i, 1,1,1 — 5*1,1,2 + 5i,i,_2 — 5i,2,i — 5i,_2,i — ^2,1,1 + 5_2,i,i) 

+^2,2 — ^2,-2 — 5_2,2 + 5-2,-2) ; (B.13) 

q(o){5)(_|) =960^5 (2(251,1,1,1,1 - 52,1,1,1 - 5i,2,i,i - 5i,i,2,i - 5i,i,i,2 + 5i,i,i,-2 - 5i,i,_2,i 

+5l, -2,1,1 — 5-2,1, 1,1) + 5i,2,2 + 52,1,2 + 52,2,1 " 5l,-2,2 — 5l,2,-2 

+5l,-2,-2 + 5-2,1,2 — 52,1,-2 — 5-2,1,-2 + 5-2,2,1 + 52,-2,1 + 5_2,-2,l) • (B.14) 

The NLO result contributes with the following terms 

g(i)"(_|) = 4i2(^a(M)(25i,i - 52 + 5_2) - 8(4(35i,i,i,i - 52,i,i - 5i,2,i - 5i,i,2) 

+2(5i,i,-2 + 5i,-2,l — 5-2, 1,1) + 5i,3 + 53,1 — 5i,_3 + 5-3,1 + 52,2 

-52,-2 + 5-2,2-5-2,-2)) (B.15) 

Q(i)'"(-§) = 24 (a(M)(25i,i,i - 5i,2 - 52,i + 5i,-2 - 5-2,i) - 8(205i,i,i,i,i 

— 7(5i,i,i,2 + 5i, 1,2,1 + 5i,2,i,i + 52, 1,1,1) + 55i,i,i,-2 + 5i, 1,-2,1 — 35i,-2,i,i 

+5-2 ,1,1,1 + 25i,2,2 + 252,1,2 + 252,2,1 — 25i,2,-2 + 5i,-2,2 — 5i,-2,-2 

— 252,1,-2 — 5-2,1,2 + 5-2,1,-2 — 5-2,2,1 — 5-2,-2,1 + 5l,l,3 + 5i,3,i + 53, 1,1 

-5i,i,-3 + 5i,-3,i - 5-3,1,1 + 53,-2 + 5-3,2)) . (B.16) 



Q^^) results in the contributions 

Q(2)(_i) = a(2)(M) + 4(65-4 + a(M)(5-2 + 52 - 25i,i)) + 2(454 - 75-3,i - 35-2,2 
-55i,-3 - 95i,3 - 52,-2 - 1052,2 - 953,1 + 2(35-2,i,i + 5i,i,-2 
+6(5i,i,2 + 5i,2,i + 52,1,1 - 45i,i,i,i))) 

Q(2)'(_i) = 2i(85-5 + a(2)(M)5i + 4a(M)(5-3 - 53 + 2(52,i + 5i,2 - 5-2,i - 35i,i,i)) 

-8(455 + 45-4,1 + 25-3,-2 + 75-3,2 - 5-2,3 - 135i,4 - 752,-3 - 2052,3 
-53,-2 - 2053,2 - 1554,1 - 45-3,1,1 + 35-2,i,2 + 35-2,2,1 + 45i,-3,i 
-45i,-2,-2 - 5i,-2,2 + 165i,i,-3 + 345i,i,3 + 75i,2,-2 + 425i,2,2 + 345i,3,i 
+752,-2,1 + 752,1,-2 — 2( — 2152,1,2 — 2152,2,1 — 1753,1,1 + 75-2,i,i,i 
+55i, -2,1,1 + 75i, 1,-2,1 + 95i,i,i,-2 + 28(5i,i,i,2 + 5i, 1,2,1 + 5i, 2,1,1 + 52,i,i,i) 

-205i,i,i,i,i))) . (B.17) 



Combining these contributions together results in ()3.27p . 
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B.2 Twist-three operators 

In analogy to the twist-two operators we choose to write the leading-order Baxter function as 

fe=0 ^ 

To obtain the expressions that contribute to the anomalous dimension we have to introduce a 
regulator, i.e. we analyze the derivatives of Q^'^\u) at u = ±z(^ + e) and take the limit e ^ 0. 
All sums now have the argument 4^. 



k=l 

M 



Q(0)'(_i)=iV^^ = 2i5i, (B.19) 



Q(0)"(_|) = i2 ^ !W ^ 4 .2 (25^^^ _ s,) , (B.20) 



k=l 



A;2 



k 

k=l 



Q(0)'"(_|) = 6 ^ - 1) = 24 (2 - - 52,i) . (B.21) 
One easily verifies that 



Q«(-§) = c(M), (B.22) 
= 2i(c(M)5i - 25i,2 - 2^2,1 - 53) . (B.23) 
According to ()B.2|1 and (|B.3p the building blocks of the anomalous dimension become 

A_(M) = -853. (B.24) 
B^{M) = 1283 + 88182, (B.25) 

such that 7^-'^) is given by ()4.8p . 

Three-loops. The NNLO anomalous dimension requires the additional derivatives of Q^^^ 



M_ 



g(0)(4)(_. ) = 24i^ ^ 52(A: - 1) =96i^( 2 (2 5i,i,i,i - ^1,1,2 - 81,2,1 - 82,1,1) 



k=i 



M_ 

' b. 



+82,2 + 83,1) , (B.26) 



Q(0){5)(_ = eoi^ ^ {8^{k - 1) - 52(fc - 1)^) = 960 (2(251,1,1,1,1 - ^1,1,1,2 

k=l 

— 81,1,2,1 — 5*1,2,1,1 — 'S'2, 1,1,1) + 5*1,2,2 + 5*2,1,2 + 5*2,2,1 + 5*1,3,1) . (B.27) 

Higher derivatives of the NLO Baxter function are given by 

Q{i)"(-i) =4i2((25i,i - 82)c{M) - 8(5i,i,2 + 5i,2,i + 52,i,i) 

+45i,3 + 10^3,1 + 8^2,2 - 354) , (B.28) 

g(i)'"(-|) = 24i3((2Si,i,i - 5i,2 - 52,i)c(M) - 12(^1,1,1,2 + 5i,i,2,i + 5i,2,i,i + ^2,1,1,1) 
+2(5S'i,i,3 + 65i,3,i + 53,1,1 + 65i,2,2 + 6^2,1,2 + 652,2,i) - 35i,4 + ^4,1 
-752,3 - 353,2) . (B.29) 
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The three-loop correction contributes with the terms 
Q(2)(-§)=c(2)(M), 

q(2)'(_|) = 2i(c(2)(M)5i - c{M){S3 + 2^1,2 + 2^2,1) + 2(ll5i,4 + 1954,i - 3^2,3 + 5^3,2) 

+8(5i,2,2 + ^2,1,2 + ^2,2,1 - 25i,3,l " 4S3,l,l) - IS^s) . (B.30) 



Plugging all terms into the expansion formula for 7^^) and transforming to the canonical basis 
one obtains ()4.13p . 

C A general method for solving an inhomogeneous Baxter 
equation 

In this appendix we formulate a method for solving any consistent inhomogeneous Baxter 
equation. For simplicity we confine ourselves to the case of twist-two operators, where the 
transfer matrix is unambiguously determined, but the generalization to more complicated 
cases (upon knowing the solution to the homogeneous equation) should be straightforward. 

Lemma 1. A minimal polynomial solution to the inhomogeneous Baxter equation of the 
form 

M-l 

{u + {fQ^'\u + i) + {u- §)' QW(n -i)- tf\u)Q^^\u) = J] afc(i ± iu)k , (C.l) 

k=0 

with the twist- two transfer matrix ^'^^(n) = 2u^ — (M^ + M + ^) and arbitrary coefficients 
Ok that are independent of the spectral parameter u is given by 

fc=0 

with 

. _ ^ (i - 1)! (-l)^r(fc + M + l) 

^'^-Z. ,2i?(,-M) ' (A:!)2r(M + l-fc) • ^^-^^ 

and Q'^^^{u) being the solution to the homogeneous equation with the coefficient in front of 
the highest power of u normalized to one. 

The proof is easily obtained using the methods of Appendix E 

This lemma allows to write down the solution to the class III terms discussed in section O 
as well as to all other inhomogeneities considered in this paper. It must be noted, however, 
that this method does not lead necessarily to the simplest representation of the solution. 

Below we will show how to bring the above-mentioned inhomogeneities of class III to 
a form, for which the Lemma 1 is directly applicable. The terms in question are a real 
combination of the following function 



+ 
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Expanding around u = allows one to check that A(u) is a polynomial of degree degA(u) = 
degQ^'^^u) — 1. Moreover the real combination A{u) + A*{u) = B{u) is a real polynomial 
of degree degi?(u) = {(iegQ^^\u) — 2). In order to find the closed form of A{u) in ()C.4|) we 
note that 

Q(0)(n)=.(. + i)M(M + l)4F3fl-^' ^! +3 J " ^ 



2, 2, 2 



1, 



(C.5) 



The key fact that makes A{u) a polynomial is that 7^'^) can be written in a similar fashion, 
since 

1-M, M + 2, 1, 1 



2, 2, 2 



(C.6) 



This allows to write A{u) as 



A{u) 



7 



(0) 



-M(M + l)4i^3 



1-M, M + 2, i + m, 1 



2, 2, 2 



iM(M + 1) ^ (1 - M)fc(M + 2)k /,i 



((l + ^^^)fc + (i-^«)fc-2(l)fc) . (C.7) 



2 A:=0 

Subsequently, using the following identities 

k k 



(C.8) 



n=0 



n=0 



where si and S2 denote Stirling numbers of, respectively, first and second kind, one finds the 
final result 



A{u) 



7 



(0) 



-M(M + 1)4^3 



1-M, M + 2, i + iu, 1 



2, 2, 2 



A/ 



1 +5^ afc(^ 



(C.9) 



fc=0 



The coefficients in the above formula are given by 



?1=/C + 1 



j = k + l 



^ ((-l)m-l _ 1^ 

E ^(-jV(-ir-'^S2(m-l,A;). 



m=fc+l 



ml 



(C.IO) 



The hypergeometric part of (|C.9p can be treated with the presented method of orthogonal 
polynomials. On the other hand, the explicit form of the coefficients in (jC.lOp allows for 
the application of lemma 1. The final result is given by (13.230 together with ()3.24p . 
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